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1 Introduction 



A distinctive feature of the topology of two-dimensional space is that a pointlike object 
renders it multiply connected. A closed path taking a particle around the object is charac- 
terized by a winding number, showing how many times the object has been encircled; and 
in the path-integral formalism, paths with different winding numbers may carry different 
overall phase factors. 

A simple physical incarnation of this general principle is a system of charged particles 
in the presence of "magnetic impurities" , i.e. magnetic flux points (two-dimensional projec- 
tions of flux tubes) at given points in space. There is no classical interaction force, since the 
magnetic field outside an impurity vanishes; however, the electromagnetic vector potential 
being topologically nontrivial, the dynamics of the particle is influenced (the Aharonov- 
Bohm effect p]). Identifying impurities with particles themselves leads to the concept of 
anyon statistics f2\ , whereby the wave function acquires a nontrivial phase factor when the 
positions of two particles are continuously interchanged. The anyon model is not exactly 
solvable for more than two particles, because the usual representation of the multiparticle 
wave function (which is now a multivalued function of its complex arguments) in terms of 
products of single-particle functions is not possible. 

Obviously, the two-dimensional model with magnetic impurities may be viewed as a 
projection of a three-dimensional system in which the impurities are infinite solenoids, and 
the longitudinal motion is free. 

Topological models of this kind are interesting because of their relevance to the frac- 
tional quantum Hall effect which is a two-dimensional phenomenon and where the 
elementary excitations are believed to be anyonic. For a particle coupled to a random 
Poissonian distribution of magnetic impurities, calculations of the density of states aver- 
aged over the positions of the impurities have been performed with some interesting 
qualitative conclusions. The model is no less fascinating, however, as an example of quite 
a generic quantum-mechanical system in multiply connected space. 

Last but not least, there is a direct connection with the problem of winding number 
distribution of random paths j^]. In fact, the partition function of one particle in the 
presence of impurities is essentially a Fourier transform of the sequence Pmim2...mM of 
probabilities that a random closed path winds rrik times around point k for all = 1, . . . , A^. 
Since it is known from past experience [Hj that the partition function of such systems 
can be computed with precision better by several orders of magnitude (with comparable 
computational time expense) by means of finding the spectrum directly rather than by 
Monte Carlo simulating the winding number distribution, one tends to believe that finding 
the levels is the right way to go in pursuit of a good numerical, and possibly analytic, 
estimate of the probability distribution in question. 

The system of one particle coupled to one impurity is trivial (standard A-B). The 
kinetic angular momentum with respect to the impurity acquires a fractional part equal 
to the magnetic flux modulo the flux quantum, and the energy levels, which depend on 
the momentum, shift accordingly. In this paper we consider the problem of one particle 
coupled to two impurities [7]. A harmonic potential is used as the most convenient long- 
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distance regulator. There appears a new dimensionless parameter, the ratio of the distance 
between the two impurities to the harmonic length. (In the thermodynamic limit, when 
the harmonic frequency vanishes, it is the ratio of that distance to the thermal wavelength 
that becomes relevant.) Searching for the wave function as a linear combination of two- 
dimensional harmonics and imposing boundary conditions at the positions of the impurities 
boils down to the standard procedure: there arises a system of linear equations for the 
coefficients, and the levels are those values of energy at which the determinant vanishes. 

We find numerically the ground state and the first two excited states and observe how 
their energies depend on the strengths of the impurities (i.e. the values of the fiuxes they 
carry) and the distance between them. At certain conditions there is a nonperturbative 
behavior at vanishing distance, whose nature is the same as in the case of the A-B ground 
state |H]: a nonsingular wave function turns into a singular one. The present calculation 
is a starting step in the evaluation of the partition function, with the subsequent inference 
of the two-point winding number distribution. 



2 Two-dimensional harmonics and boundary condi- 
tions 

The single-particle two-dimensional harmonic oscillator Hamiltonian is {h = 1) 

^_ l_f_^ }_d_ ]_d^\ !^^2 

2m ydr"^ r dr d(p'^ J 2 

Upon introducing the dimensionless length p = ^mior and energy £ = E/uj and separating 
the wave function as ?/'(p, ip) = f{p) e'^"^, the energy eigenvalue equation Hip = Eip reduces 
to the radial equation 

Its two linearly independent solutions can be expressed in terms of the confiuent hyperge- 
ometric functions 

//.£(p)=pl'le-^M(a,6,p2), (3) 
giM=P^'^e-'^U{a,b,p'), (4) 

where 

6=|/| + 1. (5) 
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M^a^h^p^) is regular at p = but in general grows exponentially at p — oo; on the 
contrary, U{a,h,p'^) falls off at infinity but is singular at the origin. It is only in the 
special case of a = — n, where n = 0, 1, 2, . . ., that the two solutions actually coincide and 
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represent a wave function that is regular both at zero and infinite p; hence, the famihar 
two-dimensional harmonic spectrum 

S = 2n+\l\ + l, (6) 

with the integer angular momentum / = 0,±1,±2,... (the wave function being single- 
valued) . 

When the particle encircles an impurity with flux $ at the origin, one can always work 
in the singular gauge where the flux has been erased but the wave function acquires an 
extra factor of exp(ia(y9), where the "impurity strength" is, by deflnition, a = $/$0) with 
the flux quantum $o = 27r/e. Thus, moving the particle around the origin brings up a 
phase factor of exp(2ia7r); for fractional a, the wave function becomes multivalued. This 
is achieved by a shift I I + a; accordingly, the spectrum becomes 

S = 2n + \l + a\ + l , (7) 

Note that the phase factor — and the spectrum — is periodic in a with period 1; i.e., it is 
only the fraction of the flux quantum that matters. 

Let us now place two impurities on the plane: one of strength ai at the point with polar 
coordinates (pojO) and another one of strength a2 at {pQ,Tr). The most general particle 
wave function inside the circle p = po is 

oo 

^_(p,^)= Y: cifiAp)e''^. (8) 

l = — OD 

Outside the same circle, adding 27r to if amounts to encircling both impurities 1 and 2, 
which should yield a phase factor of exp[2i(ai -|- a2)vr]. The wave function, therefore, must 
be 

oo 

Mp,^)= E c^,<7H«,+„„£(p)e^('+-^+'^^)^. (9) 

l=—oo 

On the circle p = pQ we must require continuity of the function il){p,ip) and its normal, 
i.e. radial derivative. There remain, however, the boundary conditions associated with 
encircling the individual impurities; those cannot be satisfled with one global relation 
between ip_{pQ,ip) and ip+iPo^f)- Rather, one has to introduce branch cuts such that 
encircling zero, one, or both impurities would involve crossing different combinations of the 
cuts, and impose phase differences between the opposite sides of each cut. Indeed, if the two 
cuts are made as shown on Fig. Q the boundary conditions are satisfled. E.g., encircling 
impurity 1 counterclockwise yields a phase factor of exp[— 2ia2vr] x exp[2i(Q;i -|- a2)vr] = 
exp[2iQ;i7r]. Thus, the boundary conditions are 

'ijj-{po,(p) = ^p+{po,(p) (0 < V? < vr) ; (10) 
^-{po,v) = e-2^"^>+(po,<^) (tt < < 27r) , (11) 



4 



Figure 1: The two branch cuts and associated phase differences. 



or, substituting the exphcit expressions (jH))^®, 

OO OO 

E fidPo) = J2 di gi+a,+a,APo) e'('+"^+-^)^ (0 < ^ < vr) ; (12) 

Z = — OO l= — OD 

OO OO 

E fiApo) e'"^ = e-^'"^'^ E di gi+a,+a,APo) e'('+"^+"^)^ {n < ^ < 2n) . (13) 

Z=— OO Z=— OO 

Exactly the same conditions have to hold for the normal derivatives dipj- (p, (f) / dp, implying 

OO OO 

E cz fieiPo) e"^ = E 9Ua,+a,APo) e^('+"^+"^)^ (0 < ^ < n) ; (14) 

Z=— OO Z=— OO 

OO OO 

E Q/(,(po)e^'^ = e-2^-- E rf^^?U+..,.(Po)e^('+"^+-)^ (vr < ^ < 27r) (15) 



Z=— OO Z=— OO 



Ifieip) = 9fi,£ip)/9p]- These equations depend on energy as a parameter, and energy 
levels are those values of £ for which all four can be satisfied simultaneously. 

3 Numerical procedure 

In order to project out single Fourier components, one should multiply Eqs. p2|) - (fT^ 
by e~'-'"^, for integer k, and integrate over (p from to 2n. In practice, however, the / 
summation will be truncated to a finite number of terms, say 2N. It is impossible to satisfy 
the boundary conditions at all values of if by adjusting the finite number of coefficients. A 
better approach, leading to greater numerical precision, is to impose those conditions at 2N 
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uniformly distributed discrete points = (2n + l)7r/ (2A^), where n = 0, 1, 2, . . . , 2N — 1 
(cf. in]). Integration over tp then gets replaced with summation over ip^- Equations (fT^ 
and (|T3|l turn into one, 

oAr f ( \ iaiTTV- sin(ai7r) + (-1)'-^= sin(a27i-) , . s 

« sm 

as do (fT^ and (fTH|) . Eliminating the c^'s yields 

5:AfcK^^)rf^ = 0, (17) 

with the matrix element 

. _ sin(ai7r) + (-1)''^ sin(a27r) , ( \ f ( \ ( ^ 

^kl[C-) — — — [Jl,£{P0) 9l+ai+a2,e\P0) - Jl,£[PO ) 9l+ai+a2,e[P0 ) 

Zily sill 

(18) 

The / summation may be over any 2N values (not even necessarily sequential): any basis 
can be chosen for expanding the wave function, but the more adequate the choice, the 
more accurate the levels found. At least for low-lying states, the right choice is to have the 
values I + ai + a2 group as closely around zero as possible. 

In the N ^ oo limit, the denominator in Eq. p8|) gets replaced with {I + ai + a2 — k)TT, 
which is precisely what one would get if Eqs. (|T ^ - ()15p were integrated over (p rather than 
summed over 

The numerical procedure is simply to scan a desired interval of energy, for given ai, 
a2, and po, to find the values that nullify the determinant of 1 1^4^/1 1. 

There is an evident fourfold symmetry in the (ai, 0^2) space. Apart from the (ai, 02) ^ 
(0^2, ai) invariance, there is P-symmetry (reversal of signs of both a's) as well as periodicity 
(oj ^— aj + 1). Taken together, these imply invariance with respect to (ai,a2) * (1 — 
a2, 1 — «i). Hence, the spectrum is symmetric with respect to both diagonals of the square 
< «i,2 < I5 and it is enough to confine oneself to any one of the four triangles formed by 
those diagonals. 

A big improvement in precision can be gained by finding one and the same level at 
different values of and then extrapolating to A^ ^ 00. In fact, comparing the results 
for different A^ is the only reliable way to estimate the error caused by the truncation. 
Convergence in A^ is better for wave functions that are less singular at impurity positions, 
since it is difficult to represent singular functions with finite Fourier series. The leading- 
term behavior of a generic wave function near an impurity of strength a (0 < a < 1), with 
the distance from the impurity r ^ 0, is r^/^^'"^^/^'. Convergence, therefore, is expected 
to be the best at a = 1/2 and the worst near the endpoints a = and a = 1. 

This is indeed what is observed. The finite- A^ data, for A^ not too small, are found to 
be quite well described by an empiric formula 

S{N)=S + cN-\ (19) 

with the exponent z/ depending rather weakly on pq, but strongly on the a's (Table 1). 
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Po = 0.75 


Po = 1-5 


8 


V 


c 


E 


V 


c 


0.1, 0.1 


1.147 


0.10 


-0.153 


1.027 


0.09 


-0.030 


0.3, 0.2 


1.280 


0.41 


-0.357 


1.043 


0.41 


-0.077 


0.5, 0.5 


1.442 


1.01 


-0.773 


1.060 


0.99 


-0.272 


0.7, 0.0 


1.135 


0.67 


-0.039 


1.023 


0.82 


-0.029 


0.9, 0.3 


1.195 


0.33 


-0.200 


1.032 


0.38 


-0.051 



Table 1: The extrapolated ground state energy and the convergence rate, from a fit to 
Eq. (Uni). The values of = 40, 80, 160, 320, 640 were used for the fitting. 

The same technique was used in Ref. |6j in the treatment of the three-anyon problem, 
the whole numerical algorithm being rather similar to the present one. The convergence 
there is much better close to the fermionic end a = 1, where the wave functions are less 
singular; a certain symmetry between a and \ — a enabled one to altogether escape the 
region < a < 1/2 with poor convergence. The trick does not work here, however: for all 
four values of (ai,a2) equivalent due to symmetry, convergence is about the same. 

4 Results and discussion 

Shown on Fig. |21 are the energies of the three lowest states as functions of the impurity 
strengths at a fixed value of distance between impurities. The values = 80, 160, 320, 
640 were used for the extrapolation by means of fitting to Eq. (fTIUl . By qualitative analogy 
with the relative motion of two anyons, introducing the impurities with positive strength 
increases the (algebraic) mean value of the particle's kinetic angular momentum. Therefore, 
the energy increases if the angular momentum for zero impurity strengths is nonnegative, 
and decreases if it is negative. (Since in this case, eigenstates of the Hamiltonian are not 
eigenstates of the angular momentum, neither, in general, are the "correct" basis states 
at «! = a2 = 0, into which the level splits as the impurities are turned on. However, the 
statement holds as pertaining to the mean value.) 

A separate issue is the dependence of the levels on the distance between impurities. 
At Po = 0, the levels are given by Eq. ((7j) with a = ai + a^- At infinite po, one has the 
single-particle oscillator spectrum (jH)). For very small a\ and 0:2, there is a clear mapping 
between the two spectra, and with the change of po, each state will slowly interpolate into 
its counterpart. As the impurity strengths increase, levels come closer to crossing each 
other. However, as long as a\ oli^ there can be only avoided crossings, since there is no 
symmetry in the system. (Reflection with respect to the x axis reverses the signs of the 
a's.) A generic case is depicted on Fig. El (a). 

True crossings become possible for a\ = 02, when there is C2 symmetry; cf. Fig. ini(b). 
These will turn into avoided crossings when that symmetry is broken in any way — say, 
resulting from a change, no matter how small, in the position or strength of one of the 
impurities. Continuity is still maintained, however: If the change is small, so will be the 
minimal distance between the levels, and their wave functions will be mutually "flipping" 
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(a) 




(c) 

Figure 2: The energies of the ground state (a) and of the first two excited states (b,c) as 
functions of ai, 0:2, at po = 0.75. 
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in the narrow region near the avoided crossing — so that the overall picture will be little 
different from a true crossing. 

There is one more peculiarity in our second example, which arises whenever ai + 0:2 is 
an integer — in this case, 1. (The spectrum is then the same at zero and infinite po-) For 
small Po, the ground state energy behaves nonperturbatively, dSo/dpol^^^Q turning infinite. 
This is due to the fact that in this and only this case, the wave function of the ground 
state (and indeed of all states with zero angular momentum) is not singular for po = 
but becomes singular when the two impurities are split up. The same reason is responsible 
for the breakdown of perturbation theory which has been long known in pertinence to the 
A-B problem: A direct attempt to treat a weak flux as a perturbation and to build up an 
expansion in a fails for the ground state, because its singular wave function, r" with a small 
a, cannot be obtained from the regular unperturbed wave function, r^, by perturbation. A 
special treatment is needed in order to enable perturbation theory to produce the correct 
answer |H]. 

States with nonzero angular momentum, on the contrary, have dSo / d po\ p^^^Q = 0, since 
their wave functions vanish at the origin fast enough. 

One final remark is in order concerning the case of big distances. At any po ^ 1, there 
will exist states whose wave functions are situated almost totally inside the po circle, as well 
as ones (high enough in the spectrum) whose wave functions are almost totally outside. 
Thus, while it is true that any given level will, at big enough po, assume its integer value 
which it has in the absence of impurities, — at any given po, no matter how big, there exist 
infinitely many levels whose energies tend to their po = values. As is commonly the 
case when there are two competing parameters, the order of limiting transitions is crucial. 
In practice, an energy cutoff will be ensured by a finite temperature; and the presence of 
the impurities will not be felt when the distance between them is much greater than the 
thermal wavelength. 

To conclude, we have construed a numerical algorithm for the solution of the problem 
of a particle coupled to two magnetic impurities in two dimensions, and found several 
low-lying states. It would be quite straightforward to generalize this to any number of 
impurities on a circle: The integrated boundary conditions would still have the form akin 
to ()lfj|l . with just the matrix elements modified. Of course, for more than two impurities, 
placing them all on one circle is a loss of generality (not so for two impurities in the 
thermodynamic limit, when the center of the harmonic potential ceases to exist). With 
several impurities not on one circle, the boundary conditions would be more complicated, 
but in principle still solvable. 

Evaluating the partition function of the system considered here and inferring the prob- 
ability distribution of winding numbers of a particle around two fixed points on a plane is 
left for future work. 
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